Solar Forcing of the Streamflow of a Continental Scale South American River 
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Solar forcing on climate has been reported in several studies although the evidence so far remains 
inconclusive. Here, we analyze the streamflow of one of the largest rivers in the world, the Parana in 
southeastern South America. For the last century, we find a strong correlation with Sunspot Number, 
in multi-decadal time scales, with larger solar activity corresponding with larger streamflow. The 
correlation coefficient is r=0.78, significant to a 99% level. In shorter timescales we find a strong 
correlation with El Nino. These results are a step toward flood prediction, which might have large 
social and economic impacts. 

PACS numbers: 92.70.Qr,93.30. Jg 
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Evidence of a solar influence on climate has been tra- 
ditionally found on records of Northern Hemisphere tem- 
perature [H, or sea surface temperature 0, usually re- 
lated to changes in solar irradiance [1, Also, a link 
between solar activity and cloud cover through the action 
of cosmic rays has been proposed Q and widely discussed 
. Recently, a number of studies took a new approach to 
the problem, looking into paleoclimatic records of atmo- 
spheric moisture. For example, evidences were found of 
the solar influence on the Asian Monsoon 0, d, @, [3, 11 1 , 
in the drought conditions in Africa [l2| and Mexico 13 1 , 
and in general in tropical precipitation regimes [l^ . The 
influence of solar activity on regional precipitations was 
also found in experiments with a global climate model 

Here, we take a different approach to the problem, 
looking into a different climatic variable, also related to 
moisture, in a different time scale: we study the stream- 
flow of the Parana River during the last 100 years (ijj . 

River streamflows are excellent climatic indicators 
since they integrate precipitations, inflltrations and 
evapotranspiration over large areas. In particular, those 
rivers with continental scale basins smooth out local va- 
riations, and can be particularly useful to study global 
forcing mechanisms. Moreover, knowledge and/or pre- 
diction of streamflow regimes is fundamental for different 
social and economic reasons, from the prediction of floods 
and droughts to planning of agricultural or hydroener- 
getic conditions. 



The Parana is one of the largest rivers in the world: 
with a basin area of over 3,100,000 km^ and a mean 
streamflow, during the last two decades of the 20^^ cen- 
tury, of 20,600 m'^/s, the Parana is the flfth river of the 
world according to drainage area and the fourth accord- 
ing to streamflow. With its origin in the southernmost 
part of the Amazon forest, it flows south collecting water 
from the countries of Brazil, Paraguay, Bolivia, Uruguay 
and Argentina, and forms one of the mightiest deltas of 
the world before its outlet in the Plata River, a few kilo- 
meters north of the City of Buenos Aires. Due to the fact 
that, unlike other rivers of similar size like the Amazon 
or the Congo, it flows through heavily populated areas, 
and that it is navigated by overseas trade ships, it has 
one of the longest streamflow data series, which covers 
the last century. 

Here we analyze the streamflow data measured at a 
gauging station located in the city of Corrientes, 900 km 
north of the outlet of the Parana. It is measured contin- 
uously from 1904, on a daily basis. The Parana's hydro- 
logical year goes from September to August, with maxi- 
mum streamflow in the (Southern Hemisphere's) summer 
months of January, February and March. We therefore 
build our yearly series integrating the flow from Septem- 
ber to August of the next year. The data are shown in 
Fig. [Tfa), together with the trend obtained with a low- 
pass Fourier fllter with a 50 years cut-off. 

It can be seen that the flow of the Parana is larger in 
the last three decades, with a mean value almost 20% 
larger than that of the flrst seventy years of the 20*^ 
century. In particular, the streamflow during the last 
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30 years has increased in the months in which the flow 
is minimum, May to December, while the flow remains 
more or less constant during the months of maximum. 
This trend has already been noticed, and was attributed 
to Amazonian deforestation which should facilitate 
water drainage. However, the same trend is also found in 
other rivers of the region like the Iguazii, whose sub-basin 
has not undergone significant changes in land use during 
the 20*^^ Century On the other hand, this trend 

can be considered as an integral part of the large-scale 
variations of the climate system [IQj . It should be noted 
that southeastern South America is one of the principal 
regions of the globe where land surface temperature has 
been increasing since 1900 j2Q] . 



ized the three quantities by substracting the mean and 
dividing by the standard deviation of each series. 




FIG. 1: (a) Parana's annual streamflow at the Corrientes 
gauging station, (b) Yearly International Sunspot Number 
(Sn)- (c) Solar irradiance reconstruction 23]. The secular 
trends, obtained with a low-pass Fourier filter with a 50 years 
cut-off, are shown as thick lines. 

As a solar activity indicator we consider the yearly 
Sunspot Number {Sn) [Sll], which is shown in Fig. [T]^b) 
together with its trend, obtained in the same way than 
that for the streamflow. Alternatively, the irradiance re- 
construction by Wang et al. [2^ can be used as a solar 
activity indicator. It differs somehow from the sunspot 
record since irradiance depends more on bright regions 
on the solar surface. 

In Fig. [2]we show the detrended time series for stream- 
flow, Sn and the irradiance reconstruction. In all cases 
we have substracted the trend shown in Fig. [1] from the 
annual data, and we have performed an 11 year running- 
mean to smooth out the solar cycle. When plotting to- 
gether different quantities, two free parameters are usu- 
ally introduced, namely, the offset and the relative scales. 
To avoid these two artificial parameters, we have normal- 




FIG. 2: The detrended time series for the Parana's streamflow 
(full line), the sunspot number (dashed line) and the irradi- 
ance reconstruction (dotted line). The detrended series were 
obtained by subtracting from each data series the correspon- 
ding trend, shown as a thick line in Fig. [T] and were smoothed 
by an 11-yr running mean to smooth out the solar cycle. The 
series were normalized by subtracting the mean and divid- 
ing by the standard deviation, to avoid introducing arbitrary 
free parameters. The Pearson's correlation coefficient of the 
streamflow data with the sunspot number is r=0.78, and with 
the irradiance reconstruction is r=0.69. 



RESULTS 

Visual agreement between the Parana's streamflow, 
the sunspot number and the irradiance reconstruction 
shown in Fig. [2] is quite remarkable. The Pearson's cor- 
relation coefficient between streamflow and SN is r=0.78, 
and between streamflow and irradiance is r=0.69. We 
performed a t-Student test to check the significance of 
these correlations, reducing the number of effective points 
to take into account the autocorrelation of the series 
and the smoothing, and we obtained a significance level 
higher than 99.99% in both cases. We also analyzed 
the significance level of the correlations with a non- 
parametric random-phase test For the SN series, we 
obtained that only a 0.4% of the random series presented 
a correlation coefficient greater than 0.78. In this way, 
we obtain a significance level of 99.6% for the correlation 
between streamfiow and SN. For the irradiance recon- 
struction, the significance obtained with this method is 
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larger than 99.99%. 

Recently, a mechanism has been proposed for the in- 
fluence of solar activity on climate, involving the modu- 
lation of galactic cosmic rays (GCR) by the interplane- 
tary magnetic field associated with the solar wind and, 
therefore, with solar activity. In this picture, GCR would 
affect cloud formation on earth, through ionization of the 
terrestrial atmosphere. Therefore, periods of higher solar 
activity, when the interplanetary magnetic field is larger, 
and therefore less GCR hit the earth, the cloud cover 
would be smaller. For this reason, it is particularly in- 
teresting to check whether there is a particularly strong 
correlation between the Parana's discharge and GCR. 

Therefore, we have also checked the correlation with 
two other solar-activity indexes. First, we considered 
the neutron count at Climax, Colorado, available since 
1953 [2J|. Since neutrons are produced when GCR hit 



the upper atmosphere, neutron count is a direct measure 
of GCR flux. Furthermore, since GCR flux in different 
parts of the world depend only on latitude, following the 
strength of the terrestrial magnetic fleld. Climax's values 
are representative of GCR flux everywhere. 

The other index we used was the aa index, which is a 
measure of the disturbance level of the Earth's magnetic 
field based on magnetometer observations of two, nearly 
antipodal, stations in Australia and England [25|, and 
it is available since 1868 [2^. It is worth pointing out 
that the aa index follows the envelope of solar activity, 
and while Sn returns to zero at each solar minimum, aa 
minima reflect the long-term level of solar activity seen 
in Fig. [Il^b). Since the Earth's magnetic fleld, which is 
affected by the solar wind, determines how much of the 
GCR flux ultimately reaches the Earth, aa can also be 
used to test the GCR-climate hypothesis. 

In both cases, we found a correlation with Parana's 
streamflow, as expected since all indexes of solar activity 
are correlated between each other. However, the correla- 
tions shown in Fig. [5] are the largest, pointing to a more 
direct correlation with solar irradiance than with GCR. 

An important point to be stressed regards the sign 
of the relationship between solar activity and river dis- 
charge reported here, which implies that wetter condi- 
tions in this area coincide with periods of higher solar ac- 
tivity. This is in agreement with paleoclimatic studies of 
the Asian monsoon 0, [23| which report an increase in 
monsoon during periods of increased solar activity. Also, 
increases in solar activity were found to be correlated 
with increased moisture over Alaska during the Holocene 
(2^ . and similar results were found in simulations of cli- 
mate during a period of reduced solar activity known as 
the Maunder Minimum [29| . 

In contrast, studies in East Africa report severe 
droughts during phases of high solar activity and in- 
creased precipitation during periods of low solar irradi- 
ation [1^. To explain these differences it has been pro- 
posed that increased solar irradiation causes more evapo- 



ration in equatorial regions, enhancing the net transport 
of moisture flux to the Indian sub-continent via monsoon 
winds 0] . A similar mechanism was found in simulations 
with a climate model, in which enhanced solar forcing 
produces greater evaporation in relative cloud-free re- 
gions in the subtropics, and the resulting moisture then 
converges into the precipitation convergence zones 15 1. 

In the American continent, droughts in the Yucatan 
Peninsula have been associated with periods of high so- 
lar activity and even proposed to explain the Mayan 
decline [13|, in contrast with the results found here. 
This is in agreement with an inverse correlation that 
was found between the southern and northern regions 
of South America, with dry periods in the South corres- 
ponding almost in phase to humid intervals in the North 
and vice versa [sot . 

The fact that solar influence is different in different 
parts of the world is of particular importance when as- 
sessing the proposed relationship between solar activity 
and climatic change, since it points out to modifications 
in circulation patterns or other mechanisms that do not 
globally affect climatic variables like moisture, but affect 
their distribution instead. 




Year 



FIG. 3: Parana's streamflow minus the multi-decadal com- 
ponent (full line) compared with the Nifiol-|-2 index (dotted 
line), which is available from 1950. Both series were normal- 
ized by subtracting the mean and dividing by the standard de- 
viation, to avoid introducing arbitrary free parameters. The 
Pearson's correlation coefficient for the whole series is r=0.65. 



Finally, in Fig [3] we show the high frequency variations 
of the Parana's streamflow, obtained by subtracting the 
11-year-running-mean from the yearly data. Also shown 
is the Nihol-|-2 index, averaged from September to Au- 
gust next year, to coincide with the hydrological year of 
the river. El Nifiol-f 2, which is available since 1950, is a 
measure of the sea surface temperatures in the Equatorial 
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Pacific Ocean, close to the South American coast (0-10° 
S, 80-90° W) [3l|. In this case, a very good accordance 
between both curves can be seen, with a correlation co- 
efficient r=0.65 and a significance larger than 99.99%. In 
particular, the large annual discharges of 1982 and 1997 
are associated with two exceptional El Nino episodes. 

This relation between the Parana's streamflow and the 
ENSO phenomenon and, in particular, the sea surface 
temperatures in a region of the tropical Pacific, was al- 
ready used in flood predictions in this basin js^]. A 
similar correlation was also found for the Indian mon- 
soon ,3S] . 

CONCLUSIONS 

Streamflow variability of the Parana river has three 
temporal components: on the secular scale, it is probably 
part of the global climatic change, which at least in this 
region of the world is related with more humid conditions; 
on the multi-decadal time scale, we found a strong cor- 
relation with solar activity, as expressed by the Sunspot 
Number, and therefore probably with solar irradiance, 
with higher activity coincident with larger discharges; on 
the yearly time-scale, the dominant correlation is with El 
Nifio. 

These correlations can be used for flood prediction: a 
regression between Parana's streamflow (S*), the sunspot 
number flltered as in Fig. [2] (Sn/), and El Nifiol-|-2 
(iVi2), gives: 

S = (63.6±59.3)S'Ar/+(3.6±0.6)xl02iVi2 + (6.5±1.3)xl0^ 

between 1949 and 1999, with a Pearson correlation coef- 
ficient r=0.66, significant to 99.99%. 

Early flood prediction, in fact, has large social and 
economic impacts: During the last flood, in 1997, 
180,000 km^ of land were covered with water, 125,000 
people have to be evacuated, and 25 people died. In all, 
the three largest fioods of the Parana during the 20^^^ 
century caused economic losses for five thousand million 
dollars. 
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